WO 2005/067391 



PCT/IL2005/000064 



DIAGNOSTIC TEST FOR PARKINSON'S DISEASE 

FIELD OF THE INVENTION 

The present invention relates to the use of molecular markers for detection, 
5 prognosis and follow up of Parkinson's disease. 

Abbreviations: ACTB: beta actin; AD: Alzheimer's disease; ALAS1: 
aminolevulinate delta synthase 1; ALDH1A1: aldehyde dehydrogenase 1 family, 
member Al; ALS: Amyotrophic Lateral Sclerosis; ARPP: 21 -cyclic AMP- 

10 regulated phosphoprotein; DA: dopamine; DEPC: diethyl pyrocarbonate; EGLN1: 
egl nine homolog 1; EIF4G1: eukaryotic translation initiation factor 4 gamma, 1; 
GAPDH: glyceraldehyde-3 -phosphate dehydrogenase; HSPA8/HSC70/HSC54: 
chaperone heat shock 70kDa protein 8; L13a: RPL13A; LGALS9: lectin, 
galactoside binding, soluble, 9; LOC56920: semaphorin sem2; LRP6: low density 

15 lipoprotein receptor-related protein 6; MAN2B1: mannosidase, alpha, class 2B, 
member 1; MPTP: N-methyl-4-phenyl-l ? 2 5 3 3 6-tetrahydropyridine; OS: Oxidative 
stress; PARVA: parvin, alpha; pc: pars compacta; PD: Parkinson's disease; 
PENK: proenkephalin; PET: positron emission tomography; PMD: post-mortem 
delay; ROS: Reactive oxygen species; SELPLG: selectin P ligand; SKP1A: S- 

20 phase kinaseassociated protein 1A; SN: Substantia nigra; SPECT: single photon 
emission tomography; SPHK1: sphingosine kinase 1; SRPK2: SFRS protein 
kinase 2; SRRM2: serine/arginine repetitive matrix 2; TMEFF1: transmembrane 
protein with EGF-like and two follistatin-like domains 1; TRIM36: tripartite 
motif-containing 36; UCHL: 1-ubiquitin C-terminal hydrolase-Ll; UPS-ubiquitin- 

25 proteasome system; UPDRS: Unified Parkinson's Disease Rating Scale; VMAT: 
vesicular monoamine member; ZSIG11: putative secreted protein ZSIG1 1. 
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BACKGROUND OF THE INVENTION 

Parkinson's disease (PD) is a progressive neurodegenerative disorder, with a 
prevalence of 1% in the population above 65 years of age, that results in 
degeneration of dopamine neurons in the substantia nigra (SN), and a consequent 
5 striatal dopamine deficiency (Bernheimer et al., 1973). The causes and mechanism 
for the degeneration of dopaminergic neurons is still elusive. Sporadic (non- 
hereditary) Parkinson's disease constitutes the most common form of the disease. 

Both hereditary and sporadic forms of PD primarily converge in impairment 
of protein handling, catabolism and in oxidative-stress damage. 
10 Potential contributing factors to the etiology of PD include ongoing selective 

oxidative stress (OS) resulting from mitochondrial dysfunction, auto-oxidation or 
enzymatic (monoamine-oxidase) oxidation of dopamine (DA) and excessive iron 
accumulation in the SN pars compacta (pc) (Riederer et al., 1989; Youdim et aL, 
1993; Gotz et al., 1994; Jenner and Olanow, 1996; Olanow and Youdim, 1996; 
15 Youdim and Riederer, 1997; Jenner, 1998). Specifically, redox-active iron has been 
observed within the melanin containing neurons that selectively die and in the rim 
of Lewy body, the morphological hallmark of PD. 

Lewy body is composed of lipids, aggregated alpha synuclein (concentrating 
in its peripheral halo) and ubiquitinated, hyperphosphorylated neurofilament 
20 proteins (Jellinger, 2003). A number of studies (Ostrerova-Golts et al., 2000; Ebadi 
et al., 2001; Turnbull et al., 2001) have shown that alpha synuclein forms toxic 
aggregates in the presence of iron and this is considered to contribute to the 
formation of Lewy body via OS. Recently, Lewy bodies with aggregated 
ubiquitinated proteins are beginning to be viewed as a defensive measure aimed at 
25 removing toxic misfolded damaged proteins, rather than a mere general site for 
precipitated misfolded proteins (Hashimoto et al., 2004; Tanaka et al., 2004). 

Misregulation of brain iron metabolism has taken center stage in 
neurodegenerative diseases since a significant number of mutated iron metabolism 
genes have now been shown to be directly involved in neurodegeneration 
30 (Felletschin et al., 2003; Youdim and Riederer, 2004). Thus, iron redox status 
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constitutes a pivotal factor contributing to the extent of protein misfolding and 
aggregation in the ageing and disease affected brain. 

Defects in ubiquitination and proteasomal protein handling are common 
features in PD and other chronic neurodegenerative diseases such as Alzheimer's 
5 disease (AD), Amyotrophic Lateral Sclerosis (ALS) and Huntington disease, and in 
ageing (Ciechanover and Brundin, 2003; Dawson and Dawson, 2003). This may, in 
turn, lead to impairment in several cellular processes linked to ubiquitination such 
as cell cycle, processing and regulation of transcription, intracellular trafficking, 
signalling pathways and degradation of normal and damaged intracellular proteins 

10 (Ciechanover and Brundin, 2003). 

Accumulating evidence points out to a crucial role for protein misfolding and 
aggregation into protein inclusions in sporadic PD, which constitutes the most 
common form of the disease. For example, a loss of 20S proteasome alpha-subunits 
(McNaught et al., 2002; McNaught et al., 2003) and reduction in the activity of the 

15 26/20S proteasome system in SNpc of sporadic PD (McNaught et al., 2003) has 
also been reported. 

Three apparently independent gene mutations in alpha synuclein, parkin and 
ubiquitin C-terminal hydrolase-Ll (UCHL-1), which are capable of impairing the 
activity of the ubiquitin-proteasome system (UPS), have been described in rare 
20 forms of hereditary PD (Dauer and Przedborski, 2003). More recently, recessive 
mutations in the DJ-1 gene were proposed to play a role in cellular response to 
oxidative stress (Bonifati et al., 2003). However, none of these genes have been 
demonstrated to be mutated in sporadic PD, which constitutes more than 95% of 
total PD cases. 

25 An approach to explore the etiology of PD is based on the study of changes 

in the up or down regulation of gene expression, which might increase the 
vulnerability of the neurons to cell death or even cause it. A number of studies have 
reported alterations in the expressions of various genes, such as decrease in 
calcium-binding protein (28kDa calbindin-D) in the SN (Iacopino and Christakos, 

30 1990) and D 3 receptor mRNA in lymphocytes (Nagai et al., 1996) from PD patients. 
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Among the neurotoxins used to induce parkinsonism in animals, N-methyl-4- 
phenyl-l,2,3,6-tetrahydropyridine (MPTP) best replicates many of the 
neurochemical and anatomical characteristics of the Parkinsonian syndrome in 
rodents, primates and other species (Dauer and Przedborski, 2003). In spite of 
5 certain limitations intrinsic to this model, namely lack of progressive nature and 
Lewy bodies (Inclusion bodies), much has been learned about the molecular events 
that lead to dopaminergic neurodegeneration. 

Using the MPTP animal model alterations in the expression of several genes 
were observed. For example, an increase of glutamate decarboxylase mRNA was 

10 observed in a subpopulation of neurons in the putamen of parkinsonian monkeys, 
which provides further evidence that striato-pallidal GABAergic neurons are 
hyperactive in MPTP-treated parkinsonian monkeys (Soghomonian and Laprade, 
1997). An increase in Bax (a cell death effector) mRNA expression in SN, with a 
concomitant increase in Bax immunoreactivity was also detected in the same model 

15 (Hassouna at al., 1996). 

More recently, the MPTP model has been employed to asses differential 
gene expression changes in the midbrain of mice (Grunblatt et al., 2001), to obtain a 
more global picture of the series of events occurring during degeneration. The 
median density microarray employed in the study contained around 1200 genes. 

20 Thus, a limited picture of the differential gene expressions that may be altered was 
observed. Nevertheless, the study demonstrated alterations in genes related to 
oxidative-stress, inflammation, nitric oxide, glutamate excitotoxicity and 
neurotrophic factors pathways, which were up or down regulated. Additional gene 
pathways, not described previously, including cell cycle and iron metabolism 

25 regulation, apoptosis, intermediary metabolism and signal transduction have also 
been observed. 

Furthermore, an increase in interleukin-1 beta (IL-lp) mRNA, in the pro- 
death genes BAD, BAX, and BID was reported in methamphetamine treated rats 
(Yamaguchi et al., 1991). Concomitantly, there were significant decreases in the 
30 anti-death genes Bcl-2 and Bcl-XL (Cadet et al., 2001; Jayanthi et al., 2001). 
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In another animal model for PD, caused by depletion of dopamine via 6- 
hydroxydopamine (6-OHDA), an increase in the glutamate decarboxylase mRNA 
expression was found in adult rats striatum as well as in neonates. A parallel 
increase in preproenkephalin and a decrease in preprodynorphin mRNA levels was 
5 observed (Laprade and Soghomonian, 1999). 

Current accepted clinical criteria for the diagnosis of PD, such as Unified 
Parkinson's Disease Rating Scale (UPDRS) (Fahn, 1987), provide high sensitivity 
for detecting parkinsonism (Brooks, 1998). However, there are no sensitive and 
specific biochemical markers that can be used to reliably diagnose clinical and 

10 especially preclinical PD for sporadic and most familial PD cases. 

One approach in the diagnosis of PD is functional imaging, which provides a 
means of discriminating typical from atypical PD, revealing characteristic patterns 
of loss of dopaminergic function. In addition, positron emission tomography (PET) 
and single photon emission tomography (SPECT) show preserved levels of striatal 

15 metabolism and dopamine receptor binding in PD, whereas levels are reduced in the 
atypical variants (Brooks, 1998). Still these tools do not give an exact diagnosis of 
PD, and often, experts in PD changed their diagnoses infrequently during the 7.6- 
year follow-up (Jankovic et al., 2000). 

Moreover, all these diagnostic methods are able to detect subjects with PD 

20 only after nearly 70% of the neurons have been degenerated, as only at this point 
symptoms appears (Zigmond et al., 1989). This, of course, makes treatment and 
maybe even rescue of the neurons nearly impossible. In view of this, it is desirable 
to diagnose the disease at an early stage and, for this reason, it is important to 
develop an early diagnostic method. \ 

25 There are already many drugs that show neuroprotective effects in vitro and 

in vivo (Drukarch and van Muiswinkel, 2001; Grunblatt et al., 2001 and 2003; 
Roller, 2002; Olanow et al., 1996 and 1998; Riederer et al., 2000 and 2002; Soto- 
Otero et al., 2002). Additionally, gene therapy was shown to be most successful in 
delaying the neurodegenerative process (Azzouz et al., 2002; Le and Frim, 2002; 

30 Luo et al., 2002; McBride and Kordower, 2002; Monville, 2002; Olanow, 2002; 
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Senior, 2002; Tenenbaum et al., 2002). The problem is that these drugs and 
methods did not show any success in patients with PD, because the beginning of 
therapy is too late and the number of surviving neurons is too small. Therefore, an 
early diagnosis may provide a better time point for the submission of therapeutic 
5 strategies that can protect against the cell death occurring in PD and to prevent the 
progress of the disease. 

Therefore, it is crucial to treat PD patients before substantial damage of 
neurons occurs, and typical symptoms appear. Thus, there is a need for an early 
diagnostic method for detection and monitoring of PD at an early stage. 

10 

SUMMARY OF THE INVENTION 

The present invention relates, in one aspect, to the use of molecular markers 
for Parkinson's disease for the purpose of diagnosis, monitoring disease progression 
and effect of treatment. 

15 In another aspect, the present invention relates to methods for diagnosis, 

prognosis and follow up of Parkinson's disease comprising the use of measurable 
biological markers, wherein said markers are one or more genes exhibiting an 
altered expression pattern, or gene products encoded thereby. 

The genes for use as molecular markers according to the invention are genes 

20 which expression is upregulated or downregulated in PD patients, as detailed 
hereinafter. 

BRIEF DESCRIPTION OF THE FIGURES 

Fig. 1 shows the relative expression levels of 137 genes differentially 

25 expressed in PD samples relative to control samples. Only genes that met the 
criteria of being altered by a factor of 1.5 relative to control and passed the 
Wilcoxon test at the significant level of p<0.05 were included. Genes are clustered 
by their relative expression levels over the 12 samples. Expression levels are color 
coded relative to the mean: green for values less than the mean and red for values 

3 0 greater than the mean. 
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Figs, 2A-B show functional cluster analysis of genes involved in biological 
processes categorized according to Gene Ontology. Pie chart showing the 
distribution of down-regulated (2A) and up-regulated (2B) genes in PD samples 
compared to control brains. The number of altered genes in the different functional 
5 groups is indicated. Each gene was assigned a single category to avoid 
overestimation of the true size of each functional group. The total number of genes 
in (2A) is 68 and in (2B) is 69. 

Figs. 3A-3B show a heat map of differentially expressed genes that changed 
in at least five out of six PD samples by a factor equal or over 1.5. Patient sample 

10 signals were compared to the average of the control sample signals (geomean). 
Down-regulated (3A) and up-regulated (3B) genes in PD samples compared to 
control brains. Expression levels are color coded relative to the mean: green for 
values less than the mean and red for values greater than the mean. 

Figs. 4A-4L show a real-time quantitative Polymerase Chain Reaction 

15 (PCR) analysis confirmation of differentially expressed genes in PD. Following 
oligonucleotide hybridization, a selected number of genes whose expression was 
altered in PD SN, were validated in three separate brain areas, SNpc, SNr and the 
cerebellum, to detect tissue specific gene alterations. For SKP1A the values were 
adjusted to be higher than one and a log Y axis was applied. Standard curves 

20 represent relative gene expression normalized to the geometric mean of four house- 
keeping genes, as described in the Materials and Methods. Data points 
corresponding to PD and control cases are indicated. ANOVA, *p<0.05 vs. control. 

Figs. 5A-5B show an immunohistology analysis confirmation of 
differentially expressed SKP1 protein in the melanized dopaminergic neurons of the 

25 SNpc from parkinsonian brains. SKP1 immunoreactivity inside the melanin- 
containing dopamine (DA) neurons are indicated by blue inclusions (see green 
arrows). SKP1 is almost absent in PD. 

30 
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DETAILED DESCRIPTION OF THE INVENTION 

The present invention relates to the use of specific molecular or biological 
markers for detection, prognosis and follow up of treatment of Parkinson's disease. 

According to the present invention, genes with altered expression pattern in 
5 PD were found, such as for example the genes ALDH1A1, ARPP-21, HSPA8, 
SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, 
PSMA5, PSMC4, HIP2, PACE4, COX6A1, PFKP, OXCT, GBE1, UQCRC2, 
LANCL1, TRIP15, PIK3CA, PLCL1, GNG5, GNAI1, VEGF, RHOB, NR4A2, 
SCL31A2, SCP2, PIGH, ARIH2, GMPR2, PP, IKBKAP, PRKACB, PTPRN2, 

10 BCAS2, IARS, PPP1R8, SEP 15, TAF9, ZFP103, WRB, TMEM4, SMARCA3, 
FMR1, PDE6D, SGCE, AUH, SLC16A7, ATP6V1E1, UGTREL1, SEC22L1, CD9, 
CDH19, DUSP1, HSA6591, ACTR3, KIF2, TUBB2, ASPA, HELOl, C3orf4, 
CBR1, XPOT, LOC51142, NY-REN-45, SETO-2, EGLN1, EIF4EBP2, LGALS9, 
LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, 

15 ZSIG11, ITGB3BP, ITGAM, COL18A1, TM4SF9, LAMB2, HS3ST2, TSTA3, 
COL5A3, PALM, MYOM1, FLNB, HMBS, KRT2A, CSK, NUDC, HYPE, GAK, 
SIAT1, CSF1R, ICSBP1, CD22, ERCC1, DNAJB5, TRAF3, MMP9, EIF4G1, 
RPL36, SRPK1, CSNK1G2, RPS6KA1, JIK, LNK, INPP5D, TCOF1, NAPG, 
SLC19A1, ITSN1, LOC51035, PMVK, C21orf2, EFEMP2, TBL1X, APRT, SPUF, 

20 GLTSCR2, ADIR, PSCD4, CBFA2T1, CUGBP1, ING4, STAT6, ZNF239, TALI, 
TAF1 1, MXD4, RDHL, LOC51 157, LRP6, MBD3, and C9orf7. 

Thus, in one embodiment, the invention relates to the use of molecular 
markers for the detection, prognosis and follow up of treatment of Parkinson's 
disease, wherein the molecular markers are one or more of the above genes or their 

25 gene products. 

It was further found by the. inventors that some of the above mentioned genes 
show an increased expression in PD patients while other genes show a decreased 
expression, as compared to healthy individuals. Tables 3 and 4 hereinafter present 
the downregulated and upregulated genes, respectively, and their characterization 
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including, for each gene, the Gene Bank number, the full gene name and its 
function. 

In accordance to the present invention, the following genes were found to be 
upregulated, namely, to exhibit an increased level of expression in PD, and will be 
5 referred to hereinafter as "the upregulated genes": EGLN1, EIF4EBP2, LGALS9, 
LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, 
ZSIG11, ITGB3BP, ITGAM, COL18A1, TM4SF9, LAMB2, HS3ST2, TSTA3, 
COL5A3, PALM, MYOM1, FLNB, HMBS, KRT2A, CSK, NUDC, HYPE, GAK, 
SIAT1, CSF1R, ICSBP1, CD22, ERCC1, DNAJB5, TRAF3, MMP9, EIF4G1, 

10 RPL36, SRPK1, CSNK1G2, RPS6KA1, JIK, LNK, INPP5D, TCOF1, NAPG, 
SLC19A1, ITSN1, LOC51035, PMVK, C21orf2, EFEMP2, TBL1X, APRT, SPUF, 
GLTSCR2, ADIR, PSCD4, CBFA2T1, CUGBP1, ING4, STAT6, ZNF239, TALI, 
TAF11, MXD4, RDHL, LOC51 157, LRP6, MBD3, and C9orf7. 

Further in accordance with the present invention, the following genes were 

15 found to be downregulated, namely, to exhibit a decreased level of expression in 
PD, and will be referred to hereinafter as "the downregulated genes": ALDH1A1, 
ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, 
PSMA3, PSMA2, PSMA5, PSMC4, HIP2, PACE4, COX6A1, PFKP, OXCT, 
GBE1, UQCRC2, LANCL1, TRIP15, PIK3CA, PLCL1, GNG5, GNAI1, VEGF, 

20 RHOB, NR4A2, SCL31A2, SCP2,PIGH, ARIH2, GMPR2, PP, IKBKAP, 
PRKACB, PTPRN2, BCAS2, IARS, PPP1R8, SEP 15, TAF9, ZFP103, WRB, 
TMEM4, SMARCA3, FMR1, PDE6D, SGCE, AUH, SLC16A7, ATP6V1E1, 
UGTREL1, SEC22L1, CD9, CDH19, DUSP1, HSA6591, ACTR3, KIF2, TUBB2, 
ASP A, HELOl, C3orf4, CBR1, XPOT, LOC51142, NY-REN-45, SETO-2. 

25 As used herein the terms "molecular markers" and "biological markers" are 

used interchangeably to include gene products that can be identified and measured 
for diagnosing PD in samples obtained from individuals. 

The term "gene products" as used herein refers to the expression product of a 
gene by conversion of the information contained in a gene into a gene product. A 
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gene product can be the direct transcript of the gene, i.e., an RNA such as mRNA, 
tRNA, or any other type of RNA, or a protein produced by translation of a mRNA. 

As used herein, the term "genes with altered expression pattern in PD" refers 
to genes that are upregulated or downregulated in PD patients in comparison, e,g., 
5 to healthy individuals. The term "upregulation" as used herein in the context of 
upregulation of a gene refers to any process which results in an increase in 
production of a gene product, e.g. RNA or protein expressed by the gene. The term 
"downregulation" as used herein in the context of downregulation of a gene refers 
to any process which results in a decrease in production of a gene product, e.g. 

10 RNA or protein expressed by the gene. 

The level of gene expression, including the level of gene upregulation or 
downregulation, can be measured by known procedures as described hereinafter. In 
general, gene upregulation or downregulation comprises any detectable increase or 
decrease in the production of a gene product. The terms upregulated and gene 

15 activation can also mean that the observed activity relative to a baseline level is a 
statistically significant difference (i.e., increase or decrease). A difference is 
typically considered to be "statistically significant" if the probability of the observed 
difference occurring by chance (the p-value) is less than some predetermined level. 
As used herein a "statistically significant difference" refers to a p-value that is 

20 O.05, preferably <0.02 for the downregulated genes (Table 3), and more preferably 
O.005 for the upregulated genes (Table 4). In accordance with the present 
invention, the increase of expression in the upregulated genes is at least 1.5-fold, 
preferably 1-3 fold (Table 4) and the decrease of expression in the downregulated 
genes is at least 0.66-fold or below, preferably 0.3-0.66 (Table 3). 

25 In one preferred embodiment of the invention, the biological markers for 

PD are one or more downregulated genes selected from the group consisting of 
ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, 
ADH5, PSMA3, PSMA2, PSMA5, PSMC4, and gene products thereof, and/or one 
or more upregulated genes selected from the group consisting of EGLN1, 
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EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, 
SPHK1, SRRM2, ZSIG1 1, and gene products thereof. 

In one embodiment, the invention relates to a method for diagnosis, 
prognosis and/or follow up of treatment of Parkinson's disease comprising the use 
of molecular markers, wherein the molecular markers are one or more of the 
upregulated and/or downregulated genes of the invention. Preferably, the molecular 
markers are one or more downregulated genes selected from the group consisting of 
ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIMS 6 5 
ADH5, PSMA3, PSMA2, PSMA5, PSMC4, and gene products thereof, and/or one 
or more upregulated genes selected from the group consisting of EGLN1, 
EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, 
SPHK1, SRRM2, ZSIG1 1, and gene products thereof. 

The methods of the invention comprise obtaining a sample from the 
individual to be tested and detecting in the sample decreased or increased levels of 
one or more gene products specific to the upregulated and/or downregulated genes 
of the invention. 

The tested individual may be an individual suspected of having PD, an 
individual exhibiting Parkinsonian-like symptoms, a de novo PD patient, a PD 
patient under treatment with a neuroprotective or another suitable PD drug, to 
follow up progress of the disease and/or to monitor effectiveness of the treatment. 
Thus, the invention provides methods for diagnosing PD at an early stage and for 
monitoring treatment of PD patients, both familial and sporadic PD patients. 

Any suitable sample can be used for detection, but preferably blood, serum 
or biopsy samples of skin obtained from the tested individual are used. 

A variety of different assays can be utilized to detect alteration in gene 
pattern expression, including methods that detect gene transcript levels and gene 
encoded protein levels. More specifically, the diagnostic and prognostic methods 
disclosed herein involve obtaining a sample from an individual and determining the 
expression level of one or more of the upregulated and/or downregulated genes in 
the sample. Usually this determined value or test value is compared against some 
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type of reference such as the expression level in control (healthy individuals, non- 
PD patients) or a baseline value. 

In accordance with the invention, also the use of the molecular markers in 
vivo is contemplated. 

5 When the gene product to be tested is RNA extracted from the sample, the 

level of gene expression, i.e the level of gene upregulation or downregulation, can 
be measured by known procedures including, but not limited to, nucleic acid probe 
arrays, Northern blots, RNase protection assays (RPA), quantitative reverse- 
transcription PGR (RT-PCR), dot blot assays and in-situ hybridization. 

10 In one preferred embodiment, nucleic acid probe arrays are used to detect 

and quantitate gene transcript (as described in the Examples below). The arrays can 
be of different types and may include probes of varying types such as, for example, 
short-length synthetic probes (20-mer or 25-mer), full length cDNA or fragments of 
gene), amplified DNA, fragments of DNA (generated by restriction enzymes, for 

15 example) and reverse transcribed DNA. The array may be a custom array, including 
probes that hybridize to particular preselected subsequences of mRNA gene 
sequences of the genes or amplification products thereof, or a generic array, 
designed to analyze mRNAs irrespective of sequence. 

In the method using probe array, nucleic acids obtained from a test sample 

20 are usually reversed transcribed into labeled cDNA, although labeled mRNA can be 
used directly. The test sample containing the labeled nucleic acids is then contacted 
with the probes of the array and, after any labeled nucleic acid related to the tested 
genes present in the sample hybridizes to the probes, the array is typically subjected 
to one or more high stringency washes to remove unbound nucleic acids and to 

25 minimize nonspecific binding to the nucleic acid probes of the arrays. Binding of 
labeled nucleic acid is detected using any of a variety of commercially available 
scanners and accompanying software programs. For example, if the nucleic acids 
from the sample are labeled with fluorescent labels, hybridization intensity can be 
determined by, for example, a scanning confocal microscope in photon counting 

30 mode. The label can provide a signal that can be amplified by enzymatic methods or 
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other labels can be used including, for example, radioisotopes, chromophores, 
magnetic particles and electron dense particles. 

Those locations on the probe array that are hybridized to labeled nucleic acid 
are detected using a reader as commercially available. For customized arrays, the 
5 hybridization pattern can then be analyzed to determine the presence and/or relative 
amounts or absolute amounts of known mRNA species in samples being analyzed. 

In another preferred embodiment, real time reverse-transcription PGR (real 
time RT-PCR) methods can be used to determine the quantity of the gene mRNA 
present in the sample (see Example below). These methods involve measuring the 

10 amount of amplification product formed during an amplification process, for 
example, by fluorogenic nuclease assays, to detect and quantitate specific 
transcripts of the genes of interest. These assays continuously measure PGR product 
accumulation using a dual-labeled fluorogenic oligonucleotide probe as in the 
approach frequently referred to in the literature simply as the TaqMan® method. 

15 The probe used in real time PCR assays is typically a short (ca. 20-25 bases) 

polynucleotide that is labeled with two different fluorescent dyes: a reporter dye at 
the 5 f terminus of the probe and a quenching dye at the 3' terminus, although the 
dyes can be attached at other locations on the probe as well. For measuring a 
specific transcript, the probe is designed to have at least substantial sequence 

20 complementarity with a probe binding site on the specific transcript. Upstream and 
downstream PCR primers that bind to regions that flank the specific transcript are 
also added to the reaction mixture for use in amplifying the nucleic acid. 

When the probe is intact, energy transfer between the two fluorophores 
occurs and the quencher quenches emission from the reporter. During the extension 

25 phase of PCR, the probe is cleaved by the 5 ? nuclease activity of a nucleic acid 
polymerase such as Taq polymerase, thereby releasing the reporter dye from the 
polynucleotide-quencher complex and resulting in an increase of reporter emission 
intensity that can be measured by an appropriate detection system. The fluorescence 
emissions created during the fluorogenic assay is measured by commercially 

30 available detectors that comprise computer software capable of recording the 
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fluorescence intensity of reporter and quencher over the course of the amplification. 
These recorded values can then be used to calculate the increase in normalized 
reporter emission intensity on a continuous basis and ultimately quantify the 
amount of the mRNA being amplified. 
5 In the assay based on dot blots and in-situ hybridization, a sample from an 

individual being tested for PD is spotted on a support, e.g., a filter, and then probed 
with labeled nucleic acid probes that specifically hybridize with nucleic acids 
derived from one or more of the upregulated or downregulated genes of the 
invention. After hybridization of the probes with the immobilized nucleic acids on 
10 the filter, unbound nucleic acids are rinsed away and the presence of hybridization 
complexes detected and quantitated on the basis of the amount of labeled probe 
bound to the filter. 

Thus, according to the invention, gene expression pattern is determined by 
extraction of RNA from biological material, preferably blood or biopsy samples of 

15 skin. The RNA is isolated rapidly by a commercially available kit. The RNA is then 
tested through hybridization to a customized GeneChip array containing the 
selected genes and relevant house-keeping genes serving for normalization, or by 
means of real time PCR for each of the selected genes. The gene expression pattern 
is determined via comparison to the expression of positive and negative control 

20 RNA (with de-novo PD and healthy subjects, respectively). The pattern of the gene 
expression obtained via one of the techniques should be similar to the pattern 
described in Table 3 or 4 in order to define the subject as PD patient. 

In another embodiment, the gene product obtained by expression of the 
upregulated or downregulated gene is a protein, that can be detected by an antibody, 

25 or a fragment thereof, able to bind to that protein. The antibody, or a fragment 
derivative, may be detectably labeled with any appropriate marker, for example, a 
radioisotope, an enzyme, a fluorescent label, a paramagnetic label, or a free radical. 

In a preferred embodiment, the invention relates to a method for diagnosing 
occurrence of Parkinson's disease in an individual exhibiting Parkinsonian-like 

30 symptoms, comprising detecting in a sample obtained from said individual a 
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decreased level of expression of one or more of the downregulated genes, more 
preferably one or more of the genes ALDH1A1, ARPP-21, HSPA8, SKP1A, 
SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, and 
PSMC4, and/or an increased level of expression of one or more of the upregulated 
5 genes, more preferably one or more of the genes EGLN1, EIF4EBP2, LGALS9, 
LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, and 
ZSIG11. An increase of the expression of one or more of the upregulated genes by 
at least about 1.5 -fold and/or a decrease of the expression of one or more of the 
downregulated genes to about 0.66 -fold or below will indicate that the individual 

1 0 suffers from PD. 

In comparison with the current methods used for diagnosis of PD, this 
method has the advantage not to use rare genetically mutations or familial history of 
the disease, but rather use general gene expression changes which occur also in 
sporadic PD. These gene expression alterations may be caused not only as a 

15 consequence of specific genetic background, but also as a consequence of 
environmental background. Therefore, this method will detect both PD patients 
carrying gene mutations and sporadic PD patients very early in the development of 
the disease, before extensive cell death and irreversible substantial damage of 
neurons occurs. 

20 When untreated, PD progresses to total disability, often accompanied by 

general deterioration of all brain functions, and may lead to an early death. 
Treated, the disorder still impairs people in varying ways. Most people respond to 
some extent to medications. The extent of symptom relief, and how long this 
control of symptoms lasts, is highly variable. 

25 Thus, in another embodiment, the invention provides a method for 

prognosticating or monitoring the treatment of a PD patient, which comprises 
detecting in a sample obtained from said patient the level of expression of one or 
more of the dowregulated or upregulated genes of the invention, whereby an 
increased level of expression of one or more of the downregulated genes, more 

30 preferably of one or more of the genes ALDH1A1, ARPP-21, HSPA8, SKP1A, 
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SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, and 
PSMC4, and/or a decreased level of expression of one or more of the upregulated 
genes, more preferably one or more of the genes EGLN1, EIF4EBP2, LGALS9, 
LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, and 
ZSIG1 1, indicate the effectiveness of the treatment. 

As shown in the Examples below, the inventors identified global major 
differentially expressed genes in the most affected brain area in PD, the substantia 
nigra pars compacta (See Tables 3 and 4). Employing Affymetrix high density 
DNA microarray, differentially expressed genes compared to aged-matched controls 
were identified. The cerebellum of PD, an unaffected brain region, served as control 
for tissue specificity. Confirmation of gene expressions was achieved by analysis 
with quantitative real-time PGR. 

It was found that the downregulated genes belong to signal transduction, 
protein degradation (e.g. ubiquitin-proteasome subunits), dopaminergic 
transmission/metabolism, ion transport, protein modification/phosphorylation and 
energy pathways/glycolysis functional classes. The upregulated genes clustered 
mainly in biological processes involving cell adhesion/cytoskeleton, extracellular 
matrix components, cell cycle, protein modification/phosphorylation, protein 
metabolism, transcription and inflammation/stress (e.g. key iron and oxygen sensor 
EGLN1). 

One major finding was the particular decreased expression of SKP1A, a 
member of the SCF (E3) ligase complex specifically in the substantia nigra (SN) of 
sporadic parkinsonian patients. SKP1A decline was accompanied by decreased 
expression in various subunits of the 26S proteasome, in energy pathways and 
signal transduction, in parallel with a prominent increase in a number of genes with 
functional activities related to cell adhesion/cytoskeleton, extracellular matrix 
components and inflammation/stress. 

SKP1 is part of the Rbx family of RING proteins (Kamura et al, 1999) 
functioning within modular multiprotein Skpl, Cullin and a substrate-recognizing 
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F-box protein (SCF). This unit allows the formation of multiple E3 complexes, 
which, in turn, are able to recognize a wide spectrum of different protein substrates. 

SCF complexes are modular: SKP1 can interact with several F-box proteins, 
which are responsible for specific target recognition, thereby providing functional 
diversity and increasing the repertoire of proteins processed by this complex. 
Humans express only one functional SKP1 isoform (Semple, 2003). Thus, the 
decrease in its expression observed, may constitute a rate limiting factor and may 
account for the accumulation of a wide spectrum of ubiquitinated protein aggregates 
in brains of PD patients such as tyrosine hydroxylase, synphilin-1, alpha- synuclein, 
phosphorylated tau (Liani et al., 2004; Meredith et al., 2004; Zhang and Goodlett, 
2004). 

Quantitative real time PGR analysis revealed that the SN expresses 
extremely low mRNA levels of SKP1A in both normal and PD subjects. For 
example, control samples display very low expression values (0.013±0.011), 
differing by more than two orders of magnitude from Rab3B (1.80±0.70), whose 
expression was not affected in SNpc of PD. A clear trend of decreased expression of 
SKPA1A in the SNpc of PD was observed. The expression of SKP1A was not 
affected in the SNr or in cerebellum of PD. 

In addition to this selective decrement in SKP1A, a concomitant decrease in 
the 20S proteasome subunits alpha-5 (PSMA5), alpha-3 (PSMA3) and alpha-2 
(PSMA2) mRNAs and in two subunits of the 19S regulatory complex of the 26S 
proteasome, the non-ATPase subunit 8 (PSMD8/Rpnl2) and the ATPase subunit 4 
(PSMC4/TBP7/Rpt3), was observed in the SN of PD brains, which may further 
contribute to dopaminergic neuronal damage in PD. 

These findings reveal novel players in the neurodegenerative scenario and 
provide potential targets for the development of novel drug compounds. 

The present invention thus further provides a method for screening for an 
agent useful for treating Parkinson's disease, which comprises identifying an agent 
that upregulates the expression of one or more of the downregulated genes, 
preferably the downregulated genes ALDH1A1, ARPP-21, HSPA8, SKP1A, 
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SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, 
PSMC4, and/or downregulates the level of expression of one or more of the 
upregulated genes, preferably EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, 
MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, ZSIG11. 

High throughput screening methods as known in the art can be used to 
identify compounds that will be useful for treatment of PD. 

The invention will now be illustrated by the following non-limiting 
examples. 

EXAMPLES 

Materials and Methods 

(i) Human brain tissue 

Post-mortem human brains were from the National German Brain-Net Center 
(Project no GA10) and the Netherlands Brain-Bank (Project no. 350). PD tissues 
were obtained from moderately to severe parkinsonism individuals based on the 
Hoehn & Yahr criteria in accordance with established guidelines (Hoehn and Yahr, 
1967). All the subjects were negative for AD pathology according to Braak & Braak 
(Braak and Braak, 1997) (for details see Table 1). The average age for PD and 
control is 76.6 and 77.8 years, respectively. The average post-mortem delay (PMD, 
time interval from death to brain-freezing at -80°C) for PD and control is 26.2 and 
19.8 hours, respectively. At autopsy, brains were dissected; one half was snap 
frozen in liquid nitrogen as slices and the other half was stored in neutral buffered 
formalin for histopathology. Stringent criteria were used in the case selection of 
human brain tissues employed in these studies. Human SN and cerebellum regions 
were obtained from control and PD donors. All procedures were in accordance with 
the with consensus criteria established by the German and Netherlands Brain-Bank 
systems and in accordance NIH Guide for the Care and Use of Laboratory Human 
Tissue and approved by the University of Wurzburg Ethics Committee (Wurzburg, 
Germany). 
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(ii) Total RNA extraction 

Total RNA was prepared with a phenol-guanidine isothiocyanate reagent 
(peqGOLD TriFast; PeQLabGmbH, Erlangen, Germany) (Lukiw and Bazan, 1997). 
RNA isolation reagents were prepared with 0.2 |iM filtered diethyl pyroearbonate 
(DEPC)-treated water (Fermentas Inc., Hanover, MD, USA) and used throughout 
the isolation procedure. Total RNA samples were spectrophotometrically scanned 
from 220 to 320 nm; the A260/A280 of total RNA was typically >1.9. In addition, 
formaldehyde agarose gel electrophoresis was conducted for quality control of total 
RNA. For all total RNAs extracted, the 28S/18S ratio was >1.5. Importantly, no 
significant differences in the spectral purity, rate of degradation, molecular size, or 
yield of SN and cerebellum total RNA between PD and control, were noted. Total 
RNAs was subjected to DNase-I digestion to get rid of genomic DNA residues and 
subsequently cleaned by the RNeasy Mini Kit (RNeasy Mini Kit; Qiagen Inc., 
Valencia, California, USA). 

(Hi) Array processing 

All experiments were performed using Affymetrix HG-FOCUS 
oligonucleotide arrays, as described at 

//www.affymetrix.com/support/technical/datasheets/human_datasheet.pd Total 
RNA from each sample was used to prepare biotinylated target RNA, with minor 
modifications from the manufacturer's recommendations 

(http://www.affymetrix.com 

Briefly, 10 jug of mRNA was used to generate first-strand cDNA by using a T7- 
linked oligo(dT) primer. After second-strand synthesis, in vitro transcription was 
performed with biotinylated UTP and CTP (Enzo Diagnostics), resulting in 
approximately 100-fold amplification of RNA. A complete description of 
procedures is available at: 
http://bioinf.picr.man.ac.ul^^ 

UK_v_2.pdf. The target cDNA generated from each sample was processed as per 
manufacturers recommendation using an Affymetrix GeneChip Instrument System 

(http://www.affymetrix.com/support/technical/manual/expression_manual 
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Briefly, spike controls were added to 15 jug fragmented cDNA before overnight 
hybridisation. Arrays were then washed and stained with streptavidin- 
phycoerythrin, before being scanned on an Affymetrix GeneChip scanner. A 
complete description of these procedures is available at: 

http://bioinf.picr.man.ac.uk/mbcf/downloads/GeneChip_Hyb_Wash_Scan_Protocol 
__v_2_web.pdf. Additionally, quality and amount of starting RNA was confirmed 
using an agarose gel. After scanning, array images were assessed by eye to confirm 
scanner alignment and the absence of significant bubbles or scratches on the chip 
surface. BioB spike controls were found to be present on all the arrays, with BioC, 
BioD and CreX also present in increasing intensity. A total of 17 gene chips were 
employed (12 for the SNpc and 5 for the cerebellum of PD patients and age- 
matched controls). 

(iv) Statistical analysis 

Genes were filtered using Mas 5 algorithm results. One expected difficulty 
when working with post mortem samples is the various degrees of degradation in 
the RNA preparations. Therefore, cRNA originating from these samples contain 
more 3 f ends than 5 f prime ends. Since the Affymetrix array probe sets are designed 
so that the 3' ends are selected whenever possible 
(http ://www. affymetrix. com/support/technical/technotes/hgu 1 3 3_design_technote .p 
df), the analysis of relatively high-3 f content samples becomes feasible. Probe sets 
detected by MAS 5 as Present (P) have a signal with a p-value lower than 0.04 
(determined from the probes contained in the probe set). This allows us to use 
present call with confidence. 

Samples that are partially degraded may hybridize to part of the probes in the 
probe sets leading to a high p-value (and detected as absent). If the probe set in at 
least four out of six samples was detected as present and all of the signals were 
greater than 20, we determined this probe set to be present. A list of 3517 probe sets 
representing genes with signals higher than 20 and detected as present (P) in all 
control samples, or with signals higher than 20 in all PD samples and detected as 
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present in four out of the six samples was generated from the 8763 probe sets 
contained on the array (supplement http://eng.sheba.co.il/genomics). 

262 probe sets differentiated between patient and control samples as 
determined by the Wilcoxon rank sum test (p-value< 0.05) (supplement 
http://eng.sheba.co.il/genomics). Probe sets were further filtered selecting those 
where the ratio of the average (geomean) signals of the PD samples and the average 
signals of the control samples exceeded 1.5 or was lower than 0.66. The 
Bonferroni's correction has not been applied because it would result in loss of most 
of the valid hits. However the use of the aforementioned cut-offs together with the 
statistical test provide a high stringency analysis and the possibility to focus on 
defined subset of genes engaged in a common biological process(s). 69 probe sets 
were up regulated and 68 were down regulated. 

(v) Real-time quantitative PCR validation of results 

In order to validate the microarray results, we performed quantitative real- 
time RT-PCR for mRNA samples from PD and controls in the SNpc, pars reticulate 
(SNr) and cerebellum. Total RNA (1-0.4 mg) from each sample was reverse 
transcribed with random hexamer and oligodT primers using iScriptTM cDNA 
Synthesis Kit (BioRad Laboratories, Hercules, CA, USA, 170-8890). The genes 
were normalized to the house-keeping genes: beta actin (ACTB), ribosomal protein 
LI 3a (RPL13A), aminolevulinate delta synthase 1 (ALAS1) and glyceraldehyde-3- 
phosphate dehydrogenase (GAPDH) (QuantiTectTM Gene Expression Assay, 
Qiagen Inc., Valencia, CA, USA, Hs ACTB Assay 241013 & Hs GAPD Assay 
241011), chosen after analysis according to the program geNorm (vs. 3.3; to be 
download from: http://medgen31.ugent.be/jvdesomp/genorm/) (Vandesompele et 
al., 2002; Schulz et al., 2004). The geNorm program determined the most stable 
house-keeping genes from a set of tested genes in a given cDNA sample panel and 
calculates a gene expression normalization factor for each tissue sample based on 
the geometric mean of a user-defined number of house keeping genes. Originally, 
we tested a total of six house-keeping genes, including ribosomal 18S and 
cyclophilin A, but the latter two were found not suitable. 
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These house-keeping genes were tested for their stability and found be the 
most accurate for normalization. Absence of DNA contamination was verified by 
amplifying the house-keeping gene ribosomal 18S, and running the reaction 
solution on agarose gel to confirm the absence of product. Minus RT samples tested 
simultaneously with experimental samples by quantitative RT-PCR with FAM 
consistently yielded no amplification below 35 cycles using the above protocol. 
Real-time PCR was performed in the iCycler iQ system (BioRad Co., Hercules, CA, 
USA) as described previously (Svaren et al., 2000; Ugozzoli et al., 2002). Briefly, 
30-100 ng of cDNA and gene specific primers and probes produced by 
QuantiTectTM Custome Assay (Qiagen Inc., Valencia, CA, USA,) (Table 2) were 
added to QuantiTect Probe PCR Master Mix (Qiagen Inc., Valencia, CA, USA, 
204343). Real-time PCR was subjected to PCR amplification (1 cycle at 95°C for 
15 min, 30-45 cycles at 94°C C for 15 s, annealing and detecting with FAM at 56°C 
for 30 s and extension at 76°C for 30 s). All PCR reactions were run in duplicate. 
The amplified transcripts were quantified using the comparative CT method 
analyzed with the BioRad iCycler iQ system program. Standard curves for each 
amplification product were generated from 10-fold dilutions of pooled cDNA 
amplicons, isolated from agarose gel using MinEluteTM Gel Extraction Kit (Qiagen 
Inc., Valencia, CA, USA), to determine primer efficiency and quantization. Data 
was analyzed with Microsoft Excel 2000 to generate raw expression values. The 
differences in gene expression in the different brain regions were compared using 
the analysis of variance (ANOVA), the StatView software program (Stat View 5.0. 
software, SAS Institute Inc. Gary, NC, USA) on a PC computer. 

vi) Immunohistochemical analysis 

Paraffin-embedded substantia nigra from parkinsonian patients and age- 
matched controls underwent serial coronal sections (6 micrometer thickness) by 
microtome. Deparaffinized sections were unmasked in digital decloacking chamber 
(Biockare Medical, Walnut Creek, CA, USA) using citrate buffer (pH 6.8), blocked 
with 10% normal goat serum for 2h at RT and incubated overnight at 4°C, with anti- 
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SKP1. Detection was obtained by appropriate biotinylated second antibody with 
streptavidin-peroxidase conjugate and S-(2-aminoethyl)-L-cysteine (AEC) as 
substrate, using Histostain-Plus kit (Zymed, South San Francisco CA, USA) 
according to manufacture instructions, Counter staining was performed with 
hematoxylin. 

Example 1. Functional classification of differentially expressed genes in PD 
patients. 

DNA microarray analysis was performed on postmortem SNpc from 6 
parkinsonian patients and 6 aged-matched controls to identify alteration in gene 
transcription associated with the disease. 

We confined our analysis to genes with at least a 1.5 fold differential 
expression and significance level of p<0.05. Of a total of 3517 valid probe sets 
analyzed (see Materials and Methods), 137 met the criteria. A heat map was 
performed depicting the relative expression levels of these 137 genes in PD and 
control samples (Fig. 1), of which 68 were transcriptionally down-regulated (Table 
3) and 69 were up-regulated (Table 4). 

Genes are clustered by their relative expression levels over the 12 samples. 
For negative control, five oligonucleotide hybridizations were performed in a non- 
related brain area, the cerebellum of both PD (2 samples) and control (3 samples), 
to assess tissue specificity of gene changes. Hierarchical clustering showed that 
both PD and normal cerebellum display a complete different pattern of expression 
(data not shown). 

Afterwards, the genes were classified into functional groups according to 
GeneOntology annotation tools (Dennis et al., 2003) 
(http://appsl.niaid.nih.gov/David/upload.asp). The Program allows finding 
regulation trends in group of genes organized according to biological processes, 
molecular function or cellular components, as defined by the Gene Ontology 
Consortium (http ://www.geneontology .org). 



23 



WO 2005/067391 



PCT/IL2005/000064 



It is recognized that a given gene may be assigned to more than one function 
or biological pathway. This may result in overestimation of the true size of each 
functional group. We overcame this by assignment to each gene a single functional 
class. As shown in Fig. 2 A, major gene downregulations were observed in the 
signal transduction, protein degradation, dopaminergic transmission and 
metabolism, ion transport, protein modification/phosphorylation and energy 
pathways/glycolysis functional classes in PD. 

With respect to differentially upregulated genes in PD, they clustered mainly 
in biological processes involving cell adhesion/cytoskeleton, extracellular matrix 
components, cell cycle, protein modification/phosphorylation, protein metabolism, 
transcription and inflammation/stress (Fig. 2B). 

Statistical analysis of gene categories, to identify genes that are the most 
over-represented or enriched, was accomplished with the EASE program (Hosack et 
al., 2003). Such groupings increase confidence in the results when the proportion of 
genes that change within a functional group is significantly greater than the fraction 
of such genes in the whole chip. 

Functional classification with "Ease score" lower than 0.05 were marked as 
over represented. This analysis revealed that the most prominent changes in the 
differentially downregulated genes occurred in proteolysis and peptidolysis (10%) 
biological process, with an EASE score of p<0.02 (Table 3, marked by an asterisk) 
and in molecular functions related to ion transporter and hydrolase activities 
(p<0.05, data not presented). Within this proteolytic group we observed decreased 
expression of both catalytic and regulatory subunits of the UPS such as 20S 
proteasome subunits alpha-5 (PSMA5), alpha-3 (PSMA3) and alpha-2 (PSMA2), 
two subunits of the 19S regulatory complex of the 26S proteasome, the non- 
ATPase subunit 8 (PSMD8/Rpnl2) and the ATPase subunit 4 
(PSMC4/TBP7/Rpt3). PSMD8 slightly exceeded the upper limit of significance, but 
a clear trend of reduction was seen in the six patients (three of them>1.5 fold of 
control). Since PSMD8 clusters together with the other UPS components, we 
considered it as functionally affected in PD. 
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On the other hand, statistical significant trends in biological processes and 
molecular functions with increased gene expression in PD were found for cell 
adhesion/cellular matrix (-15%) with an EASE score of p<0.005 (Table 4, marked 
by an asterisk) and for structural molecule and protein kinase activities (p<0.05, 
5 data not shown). 

In parallel, an exceptionally stringent analysis was done selecting those 
genes that display a 1.5 fold differential expression in at least five out of the six 
parkinsonian SN. Each patient sample signals were compared to the mean of the 
control sample signals (geomean). This evaluation is of major importance, since it 

10 allows identifying variations in gene expression within each of the post-mortem 
tissue samples that may derive from a combination of factors such as age, cause of 
death, gender, post-mortem interval, severity of the disease, etc. 

This approach restricted the total differentially expressed genes in the PD 
cases to only 20, of which 8 were decreased and 12 were increased respective to the 

15 control group (Fig. 3). 

Importantly, of the 8 downregulated genes, 3 (37.5 %) belong to DA 
neurotransmission and metabolism (Fig, 3A). These include cyclic AMP-regulated 
phosphoprotein (ARPP-21), solute carrier family 18 (vesicular monoamine member 
2, VMAT2) and aldehyde dehydrogenase 1 family, member Al (ALDH1A1). Two 

20 more genes are related to protein handling and degradation (S-phase 
kinaseassociated protein 1A (pl9A), SKP1A and the chaperone heat shock 70kDa 
protein 8, HSPA8). The last 3 genes are the SFRS protein kinase 2 (SRPK2), which 
participates in phosphorylation and protein modification processes, the tripartite 
motif-containing 36 (TRIM36) and the transmembrane protein with EGF-like and 

25 two follistatin-like domains 1 (TMEFF1), both of unknown function. 

The 12 upregulated genes (Fig. 3B) include: parvin, alpha (PARVA), lectin, 
galactoside-binding, soluble, 9 (galectin 9) (LGALS9), and selectin P ligand 
(SELPLG), belonging to the cell adhesion functional group; proenkephalin (PENK) 
and low density lipoprotein receptor-related protein (6LRP6), which are related to 

30 the cell signaling class; egl nine homolog 1 (C. elegans) (EGLN1), eukaryotic 
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translation initiation factor 4E binding protein 2 (EIF4BP2), mannosidase, alpha, 
class 2B, member 1 (MAN2B1) and sphingosine kinase 1 (SPHK1) are from the 
protein and lipid metabolism and phosphorylation categories; semaphorin sem2 
(LOC56920), involved in cell development, putative secreted protein ZSIG11 
(ZSIG11) and a serine/arginine repetitive matrix 2 (SRRM2), both of unknown 
functions. 

Example 2. Quantitative real-time PGR 

To confirm our findings of decreased expression of one of the major genes 
SKP1A, playing a key role in the UPS, we conducted real-time quantitative PCR for 
it and for two additional genes, HSC70 and VMAT2, all of them affected in at least 
five out of six parkinsonian patients by a factor < 0.67. 

We added two more samples to the experimental groups (Control 7 and 
Patient 7) in cases where the RNA used for the array hybridization was not adequate 
for the PCR. We analyzed three separate areas: the SNpc, SNr and the cerebellum. 
Cerebellum was chosen as an area not related to PD, for assessing specificity of 
gene changes. 

Real-time PCR analysis revealed that the SN expresses extremely low 
mRNA levels of SKP1A in both normal and PD subjects (Fig. 4). Although not 
statistically significant, a clear trend of decreased expression in the SNpc of PD was 
observed (Fig. 4A). The expression of SKP1A was not affected in the SNr or in 
cerebellum of PD (Figs. 4B,C). 

We confirmed significant decreased gene expression of HSP8, coding for 
HSC70/HSC54, both in the SNpc and SNr of parkinsonian patients (Figs. 4D,E), 
whereas no significant alterations were observed in cerebellum from PD and 
controls (Fig. 4F). HSC70/HSC54 exists in two isoforms, the truncated one lacking 
153 amino acids residues in the protein binding and variable domain (Tsukahara et 
al., 2000). The probe sets in the Affymetrix chip cannot discriminate between them, 
so the results from both the array hybridization and PCR analyses should represent 
a combined level of both types. 
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In addition to these novel findings, our array analysis also confirmed 
evidence for decreased expression of the vesicular DA transporter (VMAT2). Real- 
time PGR revealed a specific reduction of VMAT2 mRNA in the SNpc while the 
changes in the SNr were significantly less pronounced (Figs. 4G,H). These 
alterations were specific for the SN since no significant difference was seen in the 
cerebellum between the PD and control samples (Fig. 41). . 

To determine whether the reduced expressions found in both the arrays and 
in real-time PGR analysis was merely a reflection of a general reduction in mRNA 
levels in PD SN, or loss of synapses, we investigated the expression of three genes: 
Rab3b GTPase (Figs. 4J,K,L), syntaxin 6 and the coatomer protein complex 
subunit zeta 2 (COP zeta 2) (data not presented), which are regulators of vesicle 
trafficking from endoplasmic reticulum to Golgi, vesicle membrane docking and 
neuronal exocytosis of hormones and neurotransmitters (Gonzalez and Scheller, 
1999; Futatsumori et al., 2000; Wendler and Tooze, 2001). Their expression did not 
differ between PD and control samples in the different brain areas. 

The stringent analysis performed in the experiments above, focusing on 
selected biological meaningful pathways in PD and the substantial agreement of 
gene changes reported in literature and confirmed by real-time PGR analysis, make 
the analysis of the present results valid. 

Example 3. Detection of SKP1 protein in the melanized dopaminergic neurons 
of the SNpc from parkinsonian and healthy control brains. 

The immunohistology analysis of SKP1 protein was carried out in samples of 
melanized dopaminergic neurons of the SNpc from parkinsonian brain and healthy 
control (as described in Materials and Methods). The results summarized in Fig. 5 
show that SKP1 protein is present in normal samples but is almost absent in PD 
samples. Thus, the results obtained in Examples 1 and 2 above for mRNA 
expression of SKP1, were confirmed at the protein level by immunohistology 
analysis of SKP1 protein. 
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Summary of the results 

The results of the present invention show, for the first time, the global gene 
expression changes in the SNpc of post-mortem parkinsonian brains as compared to 
aged-matched controls. One major novel finding is the negative regulation of an 
essential component of protein catabolism, the SKP1A gene. Its decline was 
accompanied by decreased expression in various subunits of the 26S proteasome, in 
energy pathways and signal transduction, in parallel with a prominent increase in a 
number of genes with functional activities related to cell adhesion/cytoskeleton, 
extracellular matrix components and inflammation/stress. 

Reduction in SKP1A, a component of the SCF ubiquitin ligase 
multiprotein complex, in parkinsonian SNpc 

SKP1 is part of the Rbx family of RING proteins (Kamura et al., 1999) 
functioning within modular multiprotein Skpl, Cullin and a substrate-recognizing 
F-box protein (SCF). This unit allows the formation of multiple E3 complexes, 
which, in turn, are able to recognize a wide spectrum of different protein substrates. 
SCF complexes are modular: SKP1 can interact with several F-box proteins, which 
are responsible for specific target recognition, thereby providing functional 
diversity and increasing the repertoire of proteins processed by this complex. 
Humans express only one functional SKP1 isoform (Semple, 2003). Thus, the 
decrease in its expression observed herein, may constitute a rate-limiting factor and 
may account for the accumulation of a wide spectrum of ubiquitinated protein 
aggregates in brains of PD patients such as tyrosine hydroxylase, synphilin-1, a- 
synuclein, phosphorylated tau (Liani et al., 2004; Meredith et al., 2004; Zhang and 
Goodlett, 2004). 

Quantitative real time PCR analysis revealed that the SN expresses 
extremely low mRNA levels of SKP1A in both normal and PD subjects. For 
example, control samples display very low expression values (0.013±0.011), 
differing by more than two orders of magnitude from Rab3B (1.80±0.70), whose 
expression was not affected in SNpc of PD. It is also probable that the mRNA of 
SKP1A is particularly unstable and it degrades fast in post-mortem tissue. This may 
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explain the great variation in the expression levels of the samples, although a clear 
trend of decreased expression in the SNpc of PD was observed. Because of the 
special design of the Affymetrix chip, even samples that are partially degraded, 
such as in post-mortem, may in part hybridize to the probes in the probe sets leading 
to a high p-value and considered as present. Thus, the chip strategy together with 
the stringent analysis (p<0.05, decrease in five patients by a factor of at least 1.5), 
make this result confident. The low SKP1A expression levels may explain, at least 
in part, the particular susceptibility of the SNpc to different types of stresses, such 
as increased iron concentration, enzymatic (monoamine oxidase) and non- 
enzymatic (auto-oxidation) DA metabolism, abnormal protein accumulation, 
proteasomal inhibition and neurotoxin-induced OS. In this context, the PA28 
multisubunit proteasome activator (a component of the 26S proteasome) protein 
levels were shown to be very low in the SNpc of both normal and idiopatic PD 
subjects compared to other brain areas (McNaught et al., 2003), possibly 
exacerbating the already compromised DA-containing neurons to the various stress 
insults. The expression of SKP1A was not affected in the SNr or in cerebellum of 
PD. 

A recent study has demonstrated that Parkin can also function within a novel 
SCF-like complex, along with the F-box/WD repeat protein hSel-10, responsible for 
substrate recognition, and Cullinl, which however does not include SKP1 
(Staropoli et al., 2003). According to this finding, it is possible that Parkin 
associates with adaptor proteins other than hSel-10. Such complexes would likely 
display diverse substrate specificity and may explain the several targets that have 
been reported for Parkin (Dev et al., 2003). 

Impairment in proteasome subunits 

In addition to this selective decrement in SKP1A, a concomitant decrease in 
the 20S proteasome subunits alpha-5 (PSMA5), alpha-3 (PSMA3) and alpha-2 
(PSMA2) mRNAs and in two subunits of the 19S regulatory complex of the 26S 
proteasome, the non-ATPase subunit 8 (PSMD8/Rpnl2) and the ATPase subunit 4 
(PSMC4/TBP7/Rpt3), was observed in the SN of PD brains, which may further 
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contribute to dopaminergic neuronal damage in PD. The 19S complexes are located 
at one or both extremities of the 20S proteolytic core of the proteasome (Coux et al., 
1996; Voges et al., 1999) and comprise at least 18 subunits (Glickman and 
Ciechanover, 2002). They are subdivided into two subcomplexes, the "base" and the 
5 "lid" which form the portions proximal and distal to the 20S proteolytic core 
respectively, and all together form the 26S complex. The lid consists of eight 
regulatory particle non-ATPase (Rpn) subunits, including PSMD8/Rpnl2. One 
important function of the lid is to recognize multiubiquitinated proteins and other 
potential substrates of the proteasome (Glickman et aL, 1998). The base contains 

10 three non-ATPases subunits (Rpnl, Rpn2 and Rpn 10, this last common to both base 
and lid) and six putative ATPases subunits (Rptl-6), including PSMC4/TBP7. 
These interact with non-native conformations of nonubiquitinated target proteins 
(Glickman et al., 1998; Voges et al., 1999; Strickland et aL, 2000) and one of them 
S6VTBP1, was found to bind aggregated or monomeric a-synuclein (Ghee et al., 

15 2000; Snyder et al., 2003). In the same context, it has been recently reported that 
PSMC4, also called regulatory proteasomal protein S6, specifically interacts with 
both wild type and mutant synphilin-1 (Duke et aL, 2004), another presynaptic 
protein, associated with synaptic vesicles (Ribeiro et al., 2002). Synphilin-1 was 
found associated with a-synuclein and, similarly, it accumulates in the Lewy body 

20 (Wakabayashi et al., 2000). Thus, the decreased levels of PSMC4 found in the 
present invention, together with the progressive inhibition of the UPS by the slow 
but persistent aggregation of a-synuclein, may cause their accumulation in Lewy 
bodies (Jellinger, 2003; Liani et al., 2004). The observed increase in the expression 
of the glycosaminoclycan heparan-sulphate gene in PD samples, may possible 

25 contribute as well to fibrillation and aggregation of a-synuclein, as suggested 
(Cohlberg et al., 2002). Each of the six ATPases is essential and they account for 
the ATP-requirement for proteolysis and for association of the 20S and 19S 
complexes to form the 26S complex (Ghislain et al., 1993; Gordon et al., 1993). 
Therefore, the decline in the expression of PSMC4/TBP7 in brains of Parkinsonian 

30 patients may contribute to decreased levels of 26S proteasome complex, abnormal 
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accumulation of ubiquitinated proteins and reduced rates of degradation of short- 
lived proteins such as cyclins, which in turn may induce cell defects (for review see 
(Coux et al., 1996; Voges et ah, 1999)). Indeed, we observed a decline in cyclin G 
associated kinase gene. Accumulation of cyclins and cyclin-dependent kinases, has 
been reported in post-mitotic neurons undergoing apoptosis (Padmanabhan et al., 
1999; Copani et al., 2001) and in the MPTP model of PD (Grunblatt et al., 2001). It 
has been suggested that this event may represent an attempt at cell cycle re-entry 
(Verdaguer et al., 2002). Our findings are in line with a previous report of decreased 
expression of the 20S proteasome a-subunits, but not P-subunits, decreased protein 
expression levels of some 19S subunits and functional deficits in the 26/20S 
proteasome activity in the SNpc of patients with sporadic PD (McNaught et al., 
2003). In this study however, the exact nature of the different subunits affected was 
not established, as the antibodies employed in the western blot analyses recognize 
sequences that are common to the various proteasome subunits and because of 
possible protein dimerization, leading to overestimation of the molecular weight. 
Thus, the present invention provides a wider view, revealing gene changes in the 
expression of specific proteasome components. 
Heat shock protein chaperone HSC-70 

Proteotoxic insults to cells or several stress conditions can induce up 
regulation of molecular chaperones aimed at protecting cells by assisting to correct 
folding of wild type and mutated proteins. One of such is the 70 kDa heat-shock 
cognate protein (Hsc70) (Zinsmaier and Bronk, 2001), a member of the heat-shock 
protein 70 (Hsp70) family. It is suggested that Hsc70 normally mediates cellular 
processes such as protein folding, refolding, assembly, disassembly and protein 
translocation via biological membranes. It has been shown that overexpression of 
Hsp70 reduces the amount of misfolded, aggregated a-synuclein species in vivo and 
in vitro (Klucken et al., 2004), prevents the loss of dopaminergic neurons in a 
model of PD in transgenic flies that express either the wild-type a-synuclein, as 
well as the mutant forms A30P and A53T (Auluck et al., 2002) and suppresses the 
degeneration associated with polyglutamine pathogenicity (Bonini, 2002). Recently, 
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a functional polymorphism in the 5' promoter region of HSP70-1 has been reported 
in 274 PD patients, which may increase susceptibility to PD (Wu et al., 2004). We 
have found significant decreased gene expression of HSP8, coding for HSC70/ 
HSC54, both in the SNpc and SNr of five out of six parkinsonian patients, as 
confirmed by real-time quantitative PGR, whereas no significant alterations were 
observed between cerebellum from PD patients and controls, suggestive of tissue 
specificity. HSPA8 exists in two isoforms, the truncated one lacking 153 amino 
acids residues in the protein binding and variable domain. (Tsukahara et al., 2000) It 
is thought to function as an endogenous inhibitory regulator of Hsc70 by competing 
the co-chaperones (Tsukahara et al., 2000). At present, we can not differentiate 
which of these are more avidly affected. 

In addition to the involvement of Hsc70 in vesicular and nuclear trafficking 
they may play a role in protein ubiquitination through recognition of unfolded or 
aberrant proteins and delivery to a cochaperone, E3 ligase enzyme CHIP (carboxyl- 
terminus of Hsc70 interacting protein) (Murata et al., 2003). CHIP can cooperate 
with Hsp90 and/or Hsp70/ Hsc70 and ubiquitinate their attached misfolded 
substrates. Thus the Hsc70-CHIP pair represents an E3 ligase for specifically 
recognizing unfolded proteins presented by the chaperones. This molecular 
chaperone-UPS quality control system is of extreme importance in 
neurodegenerative diseases, since proteins that evade refolding or degradation by 
the UPS form aggregates that accumulates into inclusion bodies. Indeed, molecular 
chaperones as well as ubiquitin and proteasome are recruited to inclusion bodies 
and Lewy Body (Stenoien et al., 1999; Sherman and Goldberg, 2001) indicating the 
attempt of the quality control system to degrade damaged proteins or prevent their 
removal. In fact, a number of neurodegenerative diseases appear to result from 
failure of the protein quality control system (Sherman and Goldberg, 2001). 

Cell adhesion molecules, iron and oxidative stress 

An association between genes playing essential roles in neuronal 
development, such as cell migration and axonal elongation as well as components 
of the cytoskeleton and AD, has recently been suggested (De Ferrari and Inestrosa, 
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2000). A similar hypothesis can be assigned to PD where abnormal phosphorylation 
of cytoskeleton components (e.g. neurofilaments, microtubule-associated proteins) 
and ubiquitination of synaptic and other proteins, are encountered within the Lewy 
bodies (Jellinger, 2003). 

Our present findings showing abnormal upregulation of genes with structural 
and phosphorylative molecular functions integral to membrane and extracellular 
matrix (cell adhesion process), support this assumption and extend our knowledge 
revealing potential new gene targets for future therapeutics. 

During brain inflammation such as those reported in PD, AD and multiple 
sclerosis, the breakdown in the organization of key cell adhesion molecules and 
reduction in their signal transduction pathways may generate brain pathology, 
probably related to the recruitment of glia cells and macrophages and to the 
elevation in cytokines and OS. In addition, in five out of the six patients we 
observed a striking induction by more than 1.5 fold of EGLN1 (egl nine homolog 1) 
gene, a recently described proline hydroxylase enzyme belonging to the iron- and 2- 
oxoglutarate-dependent dioxygenase superfamily (Epstein et al., 2001). These 
enzymes act as key iron and oxygen sensor controlling the expression of the 
transcription factor hypoxia-inducible factor- 1 alpha (HIF), a master regulator 
orchestrating the coordinated induction of an array of hypoxia-sensitive genes. The 
target genes of HIF are especially related to angiogenesis, cell proliferation/survival 
and glucose/iron metabolism (Lee et al., 2004). Upon high oxygen levels or iron 
overload, the EGLN hydroxylases targets HIF to proteasomal degradation. 

Interestingly, the free iron-induced proteasomal-mediated degradation of iron 
regulatory protein (IRP2) involves also activation of 2-oxoglutarate-dependent 
dioxygenases and is inhibited by iron chelators (Hanson et al., 2003; Wang et al., 
2004). Thus, it is possible that IRP2 is a substrate of EGLN1, which causes post- 
translational modification, signaling it for protein degradation. Excessive 
production of EGLN1 in the SNpc may lead to a fall in IRP2 and subsequent 
decrease in transferrin receptor (TfR) mRNA and increase in ferritin levels, both 
subjected to positive and negative transcriptional regulation by IRP2, respectively 
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(Meyron-Holtz et al., 2004; Ponka, 2004) Recent studies in knock-out mice for 
IRP2, have revealed accumulation of iron in the striatum with substantial 
bradykinesia and tremor (LaVaute et al., 2001). 

Increased expression of the iron and OS sensor protein may be directly 
responsible for the observed reduction in phosphofructokinase and also the 
angiogenic factor VEGF, both regulated by the HIF proteins (Minchenko et aL, 
2003). Phosphofructokinase is the key regulatory enzyme that controls the glucose 
flux through the glycolytic pathway. Similarly, VEGF activates genes involved in 
glucose transport and metabolism via activation of the PI3K and ras pathways. The 
expression of two major players in these pathways, ras homolog gene family, 
member B and phosphoinositide-3 -kinase, catalytic, alpha polypeptide, were also 
down-regulated in parkinsonian SNpc. This finding support previous reports in 
human PD patients using positron emission tomography (PET) analysis, 
demonstrating a decrease in glucose uptake into the SN (Berding et al., 2001) and 
give a wider view of major survival pathways affected by the disease. These 
observations, together with the additional decrease in a number of energy 
pathways/glycolysis-related genes, as well as increases in iron/OS and 
inflammatory genes, is consistent with the hypothesis that mitochondrial 
dysfunction and reactive oxygen and nitrogen species contribute to the pathogenesis 
of PD. In line with this, a recent study linked a recessive mutation in the putative 
free radical sensor DJ-1 gene, with early-onset parkinsonism (Bonifati et al., 2003). 

To our surprise, apoptosis-related gene constituted the smallest functional 
class in PD samples, questioning the relevance of program cell death in the 
neurodegenerative cascade of events occurring in the disease. The role of apoptosis 
in PD is highly controversial since the evidence for it does not correlate with the 
pathological findings and the rate of neurodegeneration. We cannot exclude the 
possibility that apoptotic changes occurred at earlier stages and they decline with 
disease progression. 

Dopamine neurotransmission and metabolism 

The disruption in these neuronal processes may be exacerbated by the 
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observed reduction in the expression of cell signaling genes, which may be coupled 
to cell adhesion protein complexes and in genes belonging to the vesicle secretory 
pathway and dopaminergic neurotransmission and metabolism in the parkinsonian 
SN. In a highly stringent analysis, sought to detect crucial alterations in at least five 
out of the six PD patients, with 1.5 fold change in gene expression relative to the 
control group, we observed major decreases in the expression of 8 genes, 3 of them 
related to DA transmission, including the vesicular monoamine transporter VMAT2 
(SCL18A2). The real-time PCR confirmatory analysis revealed a specific reduction 
of VMAT2 mRNA in the SNpc while the changes in the SNr were less pronounced. 
This finding is in agreement with previous reports on post mortem SN of controls 
and PD showing a marked reduction of VMAT2 mRNA in PD which was 
associated with a marked reduction in both dopamine transporter and VMAT2 
signal per cell in the remaining pigmented neurons (Harrington et al., 1996; Brooks, 
2003). These alterations were specific for the SN as no significant difference was 
seen in the cerebellum between the PD and control samples. Two other striking 
changes were related to the reduction in DA transmission and metabolism-related 
genes ALDH1A1 and ARPP-21 coding for aldehyde dehydrogenase (ALDH) and 
cAMP regulated phosphoprotein, respectively. ARPP-21 is specifically enriched in 
DA-innervated brain regions of the basal ganglia (e.g. caudate-putamen), and in the 
substantia nigra (Ouimet et al., 1989; Tsou et al., 1993). Similar to DARPP-32 
(dopamine and adenosine 3':5'-monophosphate-regulated phosphoprotein-32K), 
ARPP-21 is activated by DA receptor Dl, and thus may represent an index of 
functional activity of D-l neurotransmission. As far as we know, this is the first 
report of decreased ARPP-21 mRNA levels in PD SN. ALDH1A1 was found to be 
expressed highly and specifically in DA cells of the SN and ventral tegmental area 
(VTA) and to be markedly reduced in SNpc dopaminergic neurons but not in those 
of the VTA of PD brains (Gaiter et al., 2003). This is in line with observation that 
striatal ALDH activity in 6-hydroxydopamine or electrical induced lesion in rats 
(Agid et al., 1973) or cats (Duncan et al., 1972), was significantly reduced. 

ALDHs are involved in the degradation of DA to aldehyde derivatives 
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(3 ,4dihy dr oxyphenylacetaldehy de and 4-hy droxy-3 -methoxyphenylacetaldehy de) 
(Mardh and Vallee, 1986), which are then metabolized to homovanilic acid and 
dihydroxyphenylacetic acid and in detoxification of aldehydes which are highly 
reactive and neurotoxic (Hjelle and Petersen, 1983). Thus, alteration in DA 
5 transmission may alter ALDH activity and/or, conversely, changes in ALDH- 
mediated metabolism may affect DA levels in nerve cell bodies and terminal fields 
in basal ganglia and the limbic system. These proteins, in conjunction with VMAT2 
may now be considered new markers for PD. 

The overall decrease in the levels of these DA neurotransmission and 

10 metabolism-related genes observed in the present invention, may seriously 
compromise neurotransmitter storage and correlate with DA neuron dysfunction. 
This may further be exacerbated by the reduced levels of SEC22L1 mRNA, a 
member of the SEC22 family of early vesicle trafficking proteins from the 
endoplasmic reticulum to (Hay et al., 1996), and by the observed reduction in many 

15 membrane carrier transporters. Impairment of cell transport and vesicle trafficking, 
carrying a wide repertoire of intracellular and membrane proteins, may well account 
for their aggregation and cell deposition into cytoplasmic inclusions (aggresomes) 
and in the Lewy body. One simple hypothesis explaining our observations is that 
loss of DA-containing neuronal bodies and neuronal synapses may account for the 

20 reduce levels of the gene products. This assumption would predict similar gene 
expression changes in a vast repertoire of vesicle traffic or synaptic-related genes. 
However, as emerged from our microarray and real-time PGR confirmations, the 
expression of many genes related to these categories such as Rab3b, syntaxin 6 and 
COP zeta 2, was not altered in PD. 

25 Protein modification/phosphorylation 

The expression of two pre-mRNA splicing related genes was inversely and 
highly affected in each of the five or six PD samples compared to control brains: the 
SRPK2 gene coding for the SFRS protein kinase 2, was down-regulated and 
SRRM2 coding for serine/arginine repetitive matrix 2, was significantly increased. 

30 SRPK2 is a kinase highly expressed in brain with a predicted sequence for specific 
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phosphorylation of arginine/serine-rich (SR) domain-containing splicing factors, 
which in turn regulate RNA splicing in brain regions (Wang et al., 1998). 
Interestingly, SRRM2 (or SRm300) 5 is a SR domain-containing splicing factor 
functioning as a co-activator within a major complex SRm 160/300, responsible for 
the processing of a subset of constitutively spliced pre-mRNAs (Blencowe et al., 
2000). 

If SRMM2 is a substrate of SRPK2, then it seems plausible that the robust 
decline in the kinase SRPK2 and the high expression of SRMM2 in PD SNpc, leads 
to abnormal hypophosphorylated SRMM2 accumulation. This in turn may affect the 
nuclear distribution of SRml60 and/or other associated SR proteins converging 
eventually in a broad impairment in splicing activity. 

Conclusion 

In spite of the vast information regarding the pathology and etiology of PD, 
it is still premature to assert what is the primary event(s) that triggers the 
development of PD. The present invention shows a reduction in the levels of 
selective components of the UPS together with a progressive misregulation of 
extracellular matrix/cytoskeleton components, concurrent with a state of OS and 
inflammation. These series of events may act independently or cooperatively during 
the course of the disease, leading eventually to the demise of dopaminergic neurons. 
Thus, subtle alterations in the kinetics of the different affected proteins may have, 
during the decades, a cumulative effect underlying the slowly progressive 
neurodegeneration of the DA-containing neurons in PD. 
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Table 1: Summary of case selection of SNpc, SNr and cerebellum regions for the 



Gene Chip array and/or quantitative real-time RT-PCR analysis. 



Case 
iNumDer 


Braak & 
i>raaK / 
Hoehn & 


Age 


Sex 


PMD (Hr) 


Cause of death 


pi 


H&Y 3/4 


73 


F 


19 


Parkinson/ renal 
insufficiency 


P2 


H&Y 5/6 


79 


F 


17 


Parkinson/ heart 
failure 


P3 


H&Y 3/4 


87 


F 


48 


n.a. 


VA 


TJjPrV Il/A 

xlOc I JV^ 


/D 


JV1 


24 


Parkinson/ 
Pneumonia 


P5 


J. IOC X .J 2 

B&BII 


/ o 




oo 
zy 


JrarKinson/ neart 
failure 


P6 


H&Y 3 


70 


1VT 

1V1 


00 

zz 


Global heart 

XdliLilC 


P7 


H&Y 5 


65 


M 


68 


± <XLi\.lLii>KJli/ VJTlUDd.1 

heart failure 


CI 




72 


M 


23 


ventricular 
fibrillation 


C2 




86 


F 


20 


renal insufficiency 


C3 




75 


M 


21 


Toxic heart failure 


C4 




85 


M 


20.5 


Left side heart 
failure 


C5 




88 


M 


48 


acute renal failure 


C6 




68 


M 


24 


n.a. 


C7 




85 


M 


25 


heart failure 



5 Post-Mortem delay, PMD; Not available, n.a. 
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Table 2: Sequences of oligonucleotide primer pairs and probes labelled with FAM 



used for real-time quantitative PGR 



Gene (GeneBank 
accession number) 


Primer pairs 
sequence (5' 

-*•) 


Primer pairs 
sequence (5' 

-3') 


Produci 
size 
(bp) 


t Cycle 
No. 


; Reaction 
Efficiency 


18s ribosomal (V01270) 


n.a. (QuantiTect Hs_RRN18S Assay] 


l 150 


36 


95 


ACTB(NM_001101) 


n.a. (QuantiTect Hs ACTB Assay) 


150 


40 


87.9 


ALAS1 (NM_000688) 


n.a. (QuantiTect Hs_AL AS 1 Assay) 


100 


40 


80.7 


GAPDH (NM_002046) 


n.a. (QuantiTect Hs_GAPD Assay) 


130 


35 


89 


RPL13A(NM_012423) 


n.a. (QuantiTect Hs_RPL13A Assay) 


100 


35 


99.5 


VMAT2 (AI269290) 


CTG TAT GTT 
CTT TGT TCT 
GGT AGA T* A 
GCT *T* GT* 
ACT TGG G 

SEQ ID NO: 1 


TTC CAG TCT 
TGC TAA CC 

SEQ ID NO: 2 


119 


45 


88.3 


HSP 70kDa (AB034951.1) 


TGC AGT TGG 
TAT TGA TCT 
TGG TAG CTT 
GGA GTG GTT 
CGG TT 

SEQ ID NO: 3 


CTC TTG TGT 
GG* GTG TT 

SEQ ID NO: 4 


107 


36 


97.5 


RAB3B (BC005035.1) 


GGGACAATG 
CAC AAG TTA 

T AAT CAA 
ACC CAA GCT 
GCT CT 

SEQ ID NO: 5 


TGG AGG AAG 
AGA GG* GT 

SEQ ID NO: 6 


107 


40 


99.6 


SKP1A(NM_006930.1) 


CAG CAG GGC 
AGA ATA AAA 
AC GAGA AGA 
ACA GAT GGA 
AA 

SEQ ID NO: 7 


GG* GAG GCA 
AAG AAA GGA 
A 

SEQ ID NO: 8 


99 


40 


99.5 


Syntaxin6 (NM_005819.1) 


GTG GTG AAA 
GGA GAG GTA 
TTG GTG GTC 
CAGTCGATT 

SEQ ID NO: 9 


CAC AGC AAC 
AAG G*G AAG 

SEQ ID NO: 10 


118 


36 


99.9 


COP zeta 2 (NM_0 16429.1) 


TTG ACT GAA 
CAG AGT GTG 

G GAT TGG 
*GG AAA TGG 
TCT GG 

SEQ ID NO: 11 | 


TGA AAT GAA 
GGC TGT GG 

SEQ ID NO: 12 


118 


40 


98.7 



* Indicates modified nucleotide 



n.a. not Available 
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Table 3: Down regulated genes in SNpc of Pakinson's Disease 



GENE BANK 


SYMBOL 


GENENAME 


FOLD OF 






* Proteolysis and Peptidolysis 




NM_006503.1 


PSMC4 


proteasome (prosome, macropain) 26S subunit, 
ATPase, 4 


0.52 


NM_002790.1 


PSMA5 


proteasome (prosome, macropain) subunit, 
alpha type, 5 


0.47 


NM_002787.1 


PSMA2 


proteasome (twosome, macronam) subunit 
alpha type, 2 


0.65 


NM_002788.1 


PSMA3 


proteasome (prosome, macropain) subunit, 
dipiio. type, j) 


0.42 


MM OOS^nQ 9 

iNIVx KJKJ '^J J jy \^ 


HTP9 


nunungiin interacting protein z 


0.54 


NM_002570.1 


PACE4 


paired basic amino acid cleaving system 4 


0.59 


NM 006930 1 


SKP1A 


0-pilai>C J^lllaoC-aboUOld-ieU. piOCem JL/\ (Jj 1 )yr\) 

jljov (limner gic i ransmission/JMetabolism 




NM 0006RQ 1 


AT TYPM A 1 


diaenyae aenyarogenase i iamiiy, member Al 


0.35 


1N1V1 UiUJUU.l 




cyclic Aivir -regulated pnospnoprotem, 2 1 kD 
Energy Pathwavs/Glvcolvsis 


0.51 


NM_004373.1 


COX6A1 


cytochrome c oxidase subunit Via polypeptide 1 


0.58 


NM_002627.1 


PFKP 


phosphofructokinase, platelet 


0.66 


NM_000436.1 


OXCT 


3-oxoacid Co A transferase 


0.57 


NM_000158.1 


GBE1 


glucan (1,4-alpha-), branching enzyme 1 
(glycogen branching enzyme, Andersen disease, 
glycogen storage disease tvt>e IV i 


0.45 


AV727381 


UQCRC2 1 


ubiquinol-cytochrome c reductase core protein 
[I 

Signal Transduction 


3.59 


NM_006055.1 ] 


LANCL1 


LanC lantibiotic synthetase component C-like 1 
^bacterial) 


3.57 


NM_004236.1 ' 


raipi5 i 


;hyroid receptor interacting protein 15 
40 


).63 
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ViJtLiNlL r5AiNJ\ 




ix JLIN ILIN AMU, 


FOLD OF 
CONTROL 


NM 006218.1 


PIK3CA 


phosphoinositide-3 -kinase, catalytic, alpha 
polypeptide 


Kf.J 


NM_006226.1 


PLCL1 


phospholipase C-like 1 


0.48 


NM 005274 1 


VJ1N vJJ 


guanine nucleotide binding protein (G protein), 
gamma 5 


U.66 


AL049933.1 


GNAIl 


guanine nucleotide binding protein (G protein), 
alpha inhibiting activity polypeptide 1 


0.64 


AF022375.1 


VEGF 


vascular endothelial growth factor 


0.49 


AI263909 


RHOB 


ras homolog gene family, member B 


0.61 


S771 54 1 


NR4A? 


nuclear receptor subfamily 4, group A, member 
2 

Transport 


0.6 


NM_00 1860.1 


SLC31A2 


solute carrier family 31 (copper transporters), 
member 2 


0.62 


NM 004731.1 


SLC16A7 


solute carrier family 16 (monocarboxylic acid 
transporters), member 7 


0.46 


BC004443 1 


ATP6V1E1 

■L A X. \J V X J—/ X 


ATPase, H+ transporting, lysosomal 31 kDa, VI 
subunit E isoform 1 




NM_005827.1 


UGTREL1 


solute carrier family 35, member Bl SEC22 
vesicle trafficking 


0.42 


AA890010 


SEC22L1 


protein-like 1 (S. cerevisiae) 


0.55 


AI269290 


SLC18A2 


solute carrier family 1 8 (vesicular monoamine), 
member 2 

Cell Adhesion 




NM_00 1769.1 


CD9 


CD9 antigen (p24) i 


D.39 


NM_02 1153.1 


CDH19 


sadherin 19, type 2 | 
Cell Cvcle 


3.5 


NM_004417.2 ] 


DUSP1 < 


dual specificity phosphatase 1 


3.57 


BC004421.1 ] 


HSA6591 s 


sine finger protein 330 
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GENE BANK 


SYMBOL 


GENENAME 


FOLD OF 

f^OTVTT? PIT 






Cell Motilitv/Cvtoskeleton 




NM 005721.2 


ACTR3 


ARP3 actin-related protein 3 homolog (yeast) 


0.44 


NM_004520.1 


KIF2 


kinesin heavy chain member 2 


0.52 


BC004188.1 


TUBB2 


tubulin, beta, 2 

Livid Biosvnthesis 


0.62 


NM_002979.1 


SCP2 


sterol carrier protein 2 


0.58 


BC004 100.1 


PIGH 


phosphatidylinositol glycan, class H 

Nucleobase. Nucleoside. Nucleotide and 
Nucleic Acid Metabolism 


0.44 


BC000422.1 


ARIH2 


ariadne homolog 2 (Drosophila) 


0.61 


IN 1VL U 1 0 J / 0. 1 


/TDD O 


guanosine monophosphate reductase 2 
Phosphate Metabolism 


0.47 


NM_02 1129.1 


PP 


pyrophosphatase (inorganic) 

Protein Modification/Phosvhorvlation 


0.56 


NM_003640.1 


IKBKAP 


inhibitor of kappa light polypeptide gene 
enhancer in B-cells, kinase complexassociated 
protein 


0.65 


NM 002731.1 


PRKACB 


protein kinase, cAMPdependent, catalytic, beta 


0.57 


NM_002847.1 


PTPRN2 


protein tyrosine phosphatase, receptor type, N 
polypeptide 2 


0.55 


NM_003138.1 


SRPK2 


SFRS protein kinase 2 

RNA Processing 


0.46 


NM_005 872.1 


BCAS2 


breast carcinoma amplified sequence 2 


0.58 


NM_0 134 17.1 


[ARS 


isoleucine-tRNA synthetase 


0.36 


NM_002713.1 : 


PPP1R8 


protein phosphatase 1, regulatory (inhibitor) 
subunit 8 

Others and Unknown Function 


D.51 


NM_004261.1 ! 


3EP15 


1 5 kDa selenoprotein ( 
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GENE BANK 


SYMBOL 


GENENAME 


FOLD OF 
CONTROL 


NM_003 187.1 


TAF9 


TAF9 RNA polymerase II, TATA box binding 
protein (TBP)-associated factor, 32kDa 


0.66 


NM_005667.1 


ZFP103 


neuroendocrine differentiation factor 


0.62 


NM_004627.1 


WRB 


tryptophan rich basic protein 


0.6 


NM_014255.1 


TMEM4 


transmembrane protein 4 


0.6 


AI760760 


SMARCA3 


SWI/SNF related, matrix associated, actin 
dependent regulator of chromatin, subfamily a, 
member 3 


0.66 


AI743037 


FMR1 


neuroendocrine differentiation factor 


0.65 


NM_002601.1 


PDE6D 


phosphodiesterase 6D, cGMPspecific, rod, delta 


0.64 


NM_003919.1 


SGCE 


sarcoglycan, epsilon 


0.63 


NM_00 1698.1 


AUH 


AU RNA binding protein/enoyl- Coenzyme A 
hydratase 


0.6 


BF439316 


TMEFF1 


transmembrane protein with EGF-like and two 
follistatin-like domains 1 


0.54 


NM_000049.1 


ASPA 


aspartoacylase (aminoacylase 2, Canavan 
disease) 


0.56 


AL1 36939.1 


HELOl 


ELOVL iamily member 5, elongation of long 
chain fatty acids (FENl/Elo2, SUR4/Elo3- like, 
yeast) 


0.6 


M30471.1 


ADH5 


alcohol dehydrogenase 5 class III), chi 
polypeptide 


0.62 


AF161522.1 


C3orf4 


chromosome 3 open reading frame 4 


0.47 


BC002511.1 


CBR1 


carbonyl reductase 1 


0.51 


AI984005 


XPOT 


exportin, tRNA (nuclear export receptor for 
tRNAs) 




NM_016139.1 


LOC51142 


phosphorylase kinase, gamma 1 (muscle) < 


3.66 


NM_016121.1 


NY-REN-45 


potassium channel tetramerisation domain 
containing 3 


3.5 


NM_022118.1 


SE70-2 


chromosome 13 open reading frame 10 ( 


166 


NM_018700.1 ' 


rRIM36 i 


tripartite motif-containing 36 ( 


3.42 
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GENE BANK 


SYMBOL 


GENENAME 


FOLD OF 
CONTROL 


AB034951.1 


HSPA8 


heat shock 70kDa protein 8 





Genes are clustered into groups by biological function. The mean fold-change of 
each gene in PD samples relative to control samples as well as the gene symbol and 
GenBank accession number are indicated. 
5 * specific categories determined to be significantly over-represented using the 
statistical clustering program EASE, are indicated 
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Table 4: Up regulated genes in SNpc of Pakinson's Disease 



GENBANK 


SYMBOL 


GENE 


FOLD 
CHANGE 






*Cell Adhesion 




NM_014288.1 


ITGB3BP 


integrin beta 3 binding protein beta3- 
endonexin) 


1 67 


rNivi_uuuojZ.z 


AAA 


integrin, alpha M (complement component 
receptor 3, alpha; also known as CD1 lb 
(pi 70), macrophage antigen alpha 

r\rj 1 vn pn ti c\ f*} 
}j\j±y jljcjj Liu.cy 


1.54 


AF018081.1 


COL18A1 


^uiidgciij v ill ? dipud 1 


1 ^1 

1.0 / 


AI741056 


SELPLG 


selectin P ligand 


2.10 


AF065389.1 


TM4SF9 


transmembrane 4 superfamily member 9 


1.60 


AK022316.1 


PARVA 


parvin, alpha 


1.80 


X79683.1 


LAMB2 


laminin, beta 2 (laminin S) 


1.76 


NM_006043.1 


HS3ST2 


heparan sulfate (glucosamine) 3- O- 

►J ULJLi. W Ll UllOltl dot 


1.84 


U58766 


TSTA3 


tissue snecifin trpmsnlfintatmn antiapn P^^R 


1 


AI984221 


COL5A3 


collagen, type V, alpha 3 

CW// Motilitv/CvtnKkplpfnrt 


1.68 


NM_002579.1 


PALM 


naralemmin 

* UrlVllUlllll 


i *jy 


NM_003803.1 


MYOM1 


myomesin 1 (skelemin) 185kDa 


2.13 


M62994.1 


FLNB 


filamin B 5 beta (actin binding protein 278) 
Development 


1.75 


NM_000 190.1 


HMBS 


hydroxymethylbilane synthase 


1.65 


NM_000423.1 


KRT2A 


keratin 2A (epidermal ichthyosis bullosa of 
Siemens) 


2.04 


NM_020 163.1 


LOC56920 


semaphorin sem2 

Cell Cvcle 


1.98 


NM_004383.1 


CSK 


c-src tyrosine kinase 


1.59 


AF241788.1 


NUDC 


nuclear distribution gene C homolog (A. 
nidulans) 


1.52 
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GENBANK 


SYMBOL 


GENE 


FOLD 

CHANGE 


NM_007076.1 


HYPE 


Huntingtin interacting protein E 


1.54 


D88435 


GAK 


cyclin G associated kinase 

Immune Response 


1.54 


AI743792 


SIAT1 


sialyltransferase 1 (betagalactoside alpha-2,6- 
sialyltransferase) colony stimulating factor 1 


1.79 


NM 005211.1 


CSF1R 


receptor, formerly McDonough feline 
sarcoma viral (v-fms) oncogene homolog 


1.57 


AI073984 


ICSBP1 


interferon consensus sequence binding 
protein 1 


1.88 


X59350 


CD22 


CD22 antigen 

Response to Stress 


1.60 


NM_00 1983.1 


ERCC1 


excision repair crosscomplementing rodent 
repair deficiency, complementation group 1 
(includes overlapping antisense sequence) 


1.61 


AK023253.1 


DNAJB5 


DnaJ (Hsp40) homolog, subfamily B, 
member 5 

Apoptosis 


1.73 


AI721219 


TRAF3 


TNF receptor-associated factor 3 

Nucleobase, Nucleoside, Nucleotide and 
Nucleic Acid Metabolism 


1 63 


NM_004349.1 


CBFA2T1 


core-binding factor, runt domain, alpha 
suDunit z, translocated to, 1, cyclm JJ-related 


1.59 


N25915 


CUGBP1 


CUG triplet repeat, RNA binding protein 1 


1.52 


AA887083 


ING4 


inhibitor of growth family, member 4 
Transcription 


1.52 


BC004973.1 


STAT6 


signal transducer and activator of 
transcription 6, interleukin-4 induced 


1.82 


NM_005674.1 


ZNF239 


zinc finger protein 239 


1.56 


NM_003 189.1 


TALI 


T-cell acute lymphocytic leukemia 1 


1.67 


AF 118094.1 


TAF11 


TAF1 1 RNA polymerase II, TATA box 
binding protein (TBP)- associated factor, 
28kDa 


1.59 
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GENBANK 


SYMBOL 


GENE 


FOLD 

CHANGE 


AK024501.1 


MXD4 


MAX dimerization protein 4 

Protein Biosynthesis 


1.82 


NM_004994.1 


MMP9 


matrix metalloproteinase 9 (gelatinase B, 
92kDa gelatinase, 92kDa type IV 
collagenase) 


2.09 


BE966878 


EIF4G1 


eukaryotic translation initiation factor 4 
gamma, 1 


1.69 


BC005057.1 


EIF4EBP2 


eukaryotic translation initiation factor 4E 
uinuing protein z 


1.75 


1N1V1__U 1 JH- 1 H . 1 




ribosomal protein L36 


1.67 






egl nine homolog 1 (C. elegans) 

Protein Modification/Phosvhorvlation 


2.37 


AW082913 


SRPKl 


SFRS protein kinase 1 


1.53 


AL530441 


CSNK1G2 


casein kinase 1, gamma 2 


1.58 


NM_002953.1 


RPS6KA1 


ribosomal protein S6 kinase, 90kDa, 
polypeptide 1 


1.67 


T T68567 1 


MAN9R1 
J vi /v j \ z, xj l 


mannosiciase, aipna, class zd, memoer 1 


O C\C\ 

2.09 


AF181985.1 


JIK 


STE20-like kinase 

Signal Transduction 


1.66 


NM_005475.1 


LNK 


lymphocyte adaptor protein 


1.51 


NM_005541.1 


INPP5D 


inositol polyphosphate-5- phosphatase, 


1.53 


NM_006211.1 


PENK 


proenkephalin 

Transport 


2.22 


AW167713 


TCOF1 


Treacher Collins-Frances chetti syndrome 1 


1 71 


JD^UU 1 007. 1 




N-ethylmaleimide-sensitive factor attachment 
protein, gamma 


1.62 


BC003068.1 


SLC19A1 


solute carrier family 19 (folate transporter), 
member 1 


1.59 


AF064243 


ITSN1 


intersectin 1 (SH3 domain protein) 


1.56 
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GENBANK 


SYMBOL 


GENE 


FOLD 
CHANGE 






Other/Unknown 




NM_015853.1 


LOC51035 


ORF 


1.57 


NM_006556.1 


PMVK 


phosphomevalonate kinase 


1.54 


NM_004928.1 


C21orf2 


chromosome 21 open reading frame 2 


1.83 


NM_016938.1 


EFEMP2 


EGF-containing fibulin-like extracellular 
matrix proiem z 


1.93 


/tJLOj J lyy 


oKJxIVIZ 


serine/ argmine repetitive matrix 2 


1.82 


AV753028 


TBL1X 


transducin (beta)~like lX-linked 


1.65 


AA927724 


APRT 


adenine phosphoribosyltransferase 


1.51 


AA679297 


SPUF 


secreted protein of unknown function 


1.70 


NM_015710.1 


GLTSCR2 


glioma tumor suppressor candidate region 
gene 2 


1.60 


NM 022371.1 


ADIR 


torsin family 3, member A 


1.52 


NM_0 15926.1 


ZSIG1 1 


putative secreted protein ZSIG1 1 


1 79 


NMJH3385.2 


PSCD4 


pleckstrin homology, Sec7 and coiled-coil 
domains 4 


2.05 


NM_005771.1 


RDHL 


dehydrogenase/reductase (SDR family) 
member 9 


1.81 


NM_0 16202.1 


LOC51157 


zinc finger protein 580 


1.87 


AF074264 


LRP6 


low density lipoprotein receptor-related 
protein 6 


2.10 


AC005943 


MBD3 


methyl-CpG binding domain protein 3 


1.66 


AL042496 


C9orf7 


chromosome 9 open reading frame 7 


1.70 


NM_009587.1 


LGALS9 


Binds galactosides. may play a role in 
thymocyte- epithelial interactions relevant 
to the biology of the thymus. The short 
isoform acts as an eosinophil 
chemoattractant 
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GENBANK 


SYMBOL 


GENE 


FOLD 
CHANGE 


XTA/T HO 1 070 1 


CDXJXr 1 

oJrrlK.1 


Catalyzes the phosphorylation of 
sphingosine to form sphingosine 1- 
phosphate (SPP) 5 a lipid mediator with 
both intra- and extracellular functions. 
Also acts on D-erythro-sphingosine and 
to a lesser extent sphinganine, but not 
other lipids, such as D,L-threo- 
dihydrosphingosine, N 5 N- 
dimethylsphingosine, diacylglycerol, 
ceramide, or phosphatidylinositol. 





* specific categories determined to be significantly over-represented using the 
statistical clustering program EASE, are indicated 
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